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Abstract 

We construct asymptotically AdSs solutions to 5-dimensional Einstein-Maxwell the- 
ory with Chern-Simons term which are dual to 4-dimensional gauge theories, including 
N = 4 SYM theory, in the presence of a constant background magnetic field B and a 
uniform electric charge density p. For the solutions corresponding to supersymmetric 
gauge theories, we find numerically that a small magnetic field causes a drastic de- 
crease in the entropy at low temperatures. The near-horizon AdS2 x i? 3 geometry of 
the purely electrically charged brane thus appears to be unstable under the addition 
of a small magnetic field. Based on this observation, we propose a formulation of the 
third law of thermodynamics (or Nernst theorem) that can be applied to black holes 
in the AdS/CFT context. 

We also find interesting behavior for smaller, non-supersymmetric, values of the 
Chern-Simons coupling k. For k = 1 we exhibit exact solutions corresponding to 
warped AdS3 black holes, and show that these can be connected to asymptotically 
AdSs spacetime. For k < 1 the entropy appears to go to a finite value at extremality, 
but the solutions still exhibit a mild singularity at strictly zero temperature. 

In addition to our numerics, we carry out a complete perturbative analysis valid to 
order B 2 , and find that this corroborates our numerical results insofar as they overlap. 



1 This work was supported in part by NSF grant PHY-07-57702. 



1 Introduction and summary of results 



The dynamics of gauge theories at finite temperature, charge density and background electro- 
magnetic fields may be studied at strong coupling via AdS/CFT dual supergravity solutions. 
The supergravity approximation is valid for large N and large 't Hooft coupling, but may be 
consistently truncated further to Einstein/Maxwell theory in the bulk when studying these 
electromagnetic effects. Temperature arises as the solutions to this Einstein/Maxwell theory 
exhibit a horizon, while charge density and background electromagnetic fields are introduced 
by imposing boundary conditions on the bulk Maxwell field. Thermodynamics and transport 
properties in gauge theories at strong coupling may then be obtained from suitable black 
hole or black brane solutions in this relatively simple Einstein/Maxwell bulk theory. 

This program has been applied extensively to 2 + 1-dimensional gauge theory, which 
is realized holographically through 3 + 1-dimensional Einstein/Maxwell theory. Its key so- 
lution is the AdS^ black brane with electric charge p, and magnetic field B. This brane 
solution is known analytically for all p and B; its spectrum of small fluctuations may be ob- 
tained systematically, and used to compute physical quantities such as electric and thermal 
conductivities [1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12]. 

While several important condensed matter problems in an external magnetic field, such 
as the Quantum Hall Effect and high T c superconductivity, are driven by 2+ 1-dimensional 
physics, it is clearly urgent to obtain results for 3 + 1-dimensional gauge theories as well. For 
instance, strong magnetic fields are created in collisions at RHIC, giving rise to observable 
effects which have been the subject of much recent discussion, e.g., [13, 14, 15, 16, 17, 18]. 

In this paper, we shall present a systematic study of the thermodynamic properties of 
3 + 1-dimensional gauge theories with finite electric charge density p in the presence of a 
constant magnetic field B. Their holographic duals should be electrically and magnetically 
charged black brane solutions to 4 + 1-dimensional Einstein/Maxwell theory with a Chern- 
Simons term. The Chern-Simons coupling k captures the strength of the anomaly of the 
boundary current, and is required to take a specific value k = 2/a/3 (in our conventions) if 
the Einstein/Maxwell theory is to be the bosonic truncation of minimal D = 5 supergravity; 
see e.g. [19]. 2 It has been proven that any M = 1 superconformal theory with an AdSs 
supergravity dual obtained by compactification from IIB or M-theory admits a consistent 
truncation to D = 5 minimal gauged supergravity [20, 21, 22]. Thus the results we find 
pertain to a large class of theories, of which M = 4 super Yang-Mills is but one example. 

The purely electric solution (B = 0) is the Reissner-Nordstrom black hole in AdS$; 
its analytic form and thermodynamic properties are well-known. The existence of purely 
magnetic solutions (p = 0) was demonstrated numerically in a previous paper [23], but no 
analytical solutions are available at present. Purely magnetic solutions interpolate between 

2 Whenever B = 0, the Chern-Simons coupling is immaterial, and does not enter into the physical quan- 
tities considered here. 
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AdSs and AdS% x R 2 near the horizon. As a function of temperature T, their entropy density 
behaves as T 3 for large T, and vanishes as BT for small T. These limits agree with M = 4 
SYM calculations at zero gauge coupling, up to factors of 3/4 and ^4/3 respectively. On the 
M = 4 SYM side the low temperature thermodynamics is governed by fermions in the lowest 
Landau level; an appealing feature of these supergravity solutions is that they reproduce this 
low temperature behavior. 

Here we wish to extend these results to nonzero p and B. Besides their clear usefulness 
for applications of AdS/CFT, this investigation has conceptual implications for the status 
of extremal black holes, as we now pause to discuss. 



1.1 Extremal black holes, Nernst's Theorem, the third law of ther- 
modynamics, and all that 

A striking feature of the Reissner-Nordstrom black brane solution, in any dimension, is that 
it possesses a smooth zero temperature limit with nonzero entropy density. This extremal 
solution exhibits a near-horizon AdS2 x R D ~ 2 region, the existence of which has played a 
central role in recent holographic descriptions of non- Fermi liquids [24, 25, 26]. However, 
a much discussed cause for concern is that while the extremal entropy apparently plays a 
crucial role in this analysis, it is not expected from the point of view of interacting fermions, 
nor from the point of view of the dual field theory where the existence of massless charged 
bosons suggests that Bose condensation should rule. One possibility is that the large ground 
state degeneracy should be understood as an artifact of the large N limit, as discussed in [9]. 
Another is that one should focus on alternative bulk theories where the extremal entropy 
vanishes, as discussed recently in the case of gauge fields coupled to massless scalars in [10]. 
While this may be the case, we would like to propose another resolution, based on the results 
we find for the response to magnetic fields. 

The tension between the existence of extremal black hole entropy and the thermody- 
namic behavior of typical systems has been discussed periodically over the years (see e.g. 
[27]), and can be phrased in terms of a clash with Nernst's "theorem" and the third "law" of 
thermodynamics. These statements can be expressed in various ways, but essentially they 
stipulate that the entropy density s should go to zero at zero temperature (see e.g. [28]). 
Despite the name, this "law" / "theorem" is actually meant to be a phenomenological obser- 
vation characterizing the behavior of observed physical systems. Indeed, only a moment's 
thought is required to realize that it is trivial to concoct theoretical counterexamples based 
on free systems. However, such counterexamples are to be thought of as being fine-tuned to 
an unphysical degree, as any realistic system will exhibit some degree of interactions, and 
these will typically lift the ground state degeneracy. The relevant question is whether the 
ground state degeneracy is stable under adding generic weak interactions or perturbations 
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of the system. If Aj represent some set of coupling constants, we should consider 

limlim S (A,,T). (1.1) 

If this limit gives zero for "typical" couplings Aj, then we may conclude that a fine tuning is 
required to sustain the entropy. 

Phrased in this way, Nernst's "theorem" admits a natural formulation in terms of black 
holes in the AdS/CFT correspondence. We should ask whether the extremal black branes 
exhibiting finite entropy at extremality are fine-tuned in the same sense when we change 
the interactions. In AdS/CFT we change the Lagrangian of the CFT by changing boundary 
conditions for fields in the bulk, and thus we can ask whether the extremal entropy persists 
even in the presence of nontrivial boundary conditions for a suitable set of bulk fields. In a 
scenario in which the entropy is to be regarded as requiring fine-tuning, we would expect to 
see behavior like that in Fig. 1. 



s 




Figure 1: Illustration of fine tuning required to maintain nonzero extremal entropy. Small 
nonzero couplings Aj lead to no appreciable effect at high temperature, but cause the low 
temperature entropy to flow to zero. 

In this work we study the case in which Aj corresponds to turning on a constant external 
magnetic field coupling to the R-current of Af = 4 SYM (or other superconformal field 
theories described holographically) , and we will present evidence that the extremal entropy 
indeed is unstable in the above sense under the inclusion of a magnetic field. Since we 
proceed numerically, and our numerics break down at very low temperature, we are not able 
to follow the entropy all the way down to zero, but the simplest extrapolation suggests a 
picture in concordance with Fig. 1. In further support of this interpretation, we will see 
that an attempt to construct a finite entropy solution perturbatively in B breaks down at 
very low temperature. These results suggest that conclusions drawn from the extremal black 
brane solution should be viewed with caution, as they can be drastically affected by even a 
small (perhaps even arbitrarily small) magnetic field. It is an interesting question to explore 
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the effect of introducing boundary conditions for other fields and to study their effect on 
extremal branes in various dimensions [29]. In this regard, it is worth noting that although 
the asymptotically AdS4 extremal brane solution maintains its entropy in the presence of a 
magnetic field, the question remains regarding more general perturbations [29]. 

1.2 Summary of results 

One of the main results of the present paper is that the low temperature thermodynamics of 
solutions carrying nonzero charge density p and magnetic field B depends crucially on the 
value of the Chern-Simons coupling k. There are three qualitatively distinct cases: k < 1, 
k — 1, and k > 1. In our conventions, the supersymmetric value is k = 2/a/3, and so falls 
into the k > 1 category. At high temperatures there is no significant distinction between the 
three cases. However, as we take the temperature to zero, holding fixed the dimensionless 
ratio B 3 /p 2 , we find markedly different behavior, as shown in Fig. 2. 



AdS 2 x r 3 k < 1 k = 1 k > 1 




Figure 2: Schematic illustration of flows in parameter space for three ranges of k. Blue lines 
are flows at various fixed values of B 3 /p 2 . Arrows indicate direction of decreasing temper- 
ature. Red lines indicate the boundary of allowed (b, q) values where nonsingular solutions 
are possible. The near-horizon geometries at the end points of the flows are indicated in the 
k < 1 and k = 1 cases. For k > 1 there exists an AdS3 x R 2 solution at b = \^3, indicated 
by the dot, but the flows do not reach this point. In the k < 1 and k > 1 cases, near the 
endpoint of the flow B 3 /p 2 becomes a very sensitive function of (b,q), depending on the 
precise direction of approach. 

We parameterize our solutions by b and q, which represent the values of the magnetic 
field and charge density at the horizon in a particular coordinate system. They differ from 
the physical magnetic field and charge density, which we are calling B and p; the latter 
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are measured at spatial infinity, and the relation between the two sets of parameters is 
determined numerically. For nonsingular solutions, b and q take values in a bounded region, 
which we can scan over numerically. As we lower the temperature holding B 3 /p 2 fixed 3 we 
flow along lines in the (b, q) plane, in what are essentially renormalization group flows. The 
flows for the various values of k are shown schematically in Fig. 2. The figures representing 
our actual numerical data will be presented in section 7, specifically in Fig. 6. 

The flows are driven towards three distinct endpoints, depending on the value of k. 

• For k < 1 we are driven towards (b — 0, q — This is the Reissner-Nordstrom 
black brane with near horizon geometry AdS 2 x R 3 . 

• For k — 1 we flow towards the curve q 2 + 2b 2 = 6; the solutions along this curve have 
a near-horizon geometry corresponding to warped AdS 3 xR 2 . These warped AdS 3 
geometries have attracted attention recently in the context of topologically massive 
gravity [30] and the Kerr/CFT correspondence [31]; here we find that they emerge as 
solutions of 4+1 dimensional Einstein-Maxwell theory, and can be connected to asymp- 
totically AdS 5 spacetimes. As we move along the curve the near-horizon geometries 
interpolate continuously between AdS 2 x R 3 and AdS 3 x R 2 . 

• For k > 1, including the supersymmetric value, we flow towards (b c ,q = 0), where b c 
is a k dependent number that starts out at for k — 1 and then decreases with k. 
To the right of this end point, at (y/3, 0), is the AdS 3 x R 2 solution that was studied 
in [23]. At the supersymmetric value of the Chern-Simons coupling, k = 2/y/3, the 
values of q and b are bounded by a critical curve, which to about 0.5% accuracy, is 
given by the relation q 2 + ab 2 = 6 with a ~ 2.44149 (this value of a is chosen to give 
high precision at the endpoint of the curve.) 

As would be expected from the flow diagrams, the behavior of the entropy at low temper- 
ature depends on k. In Fig. 3 we show our numerical results for the supersymmetric value 
k = 2/y / 3 with B 3 / p 2 fixed at and approximately .15. We plot dimensionless versions 
of the entropy density and temperature, since the dimensionful versions have no intrinsic 
meaning. 

This plot illustrates that a small value of B causes a large decrease in the entropy. This 
behavior is representative of the k > 1 case in general, and the effect seems to get more 
pronounced with increasing k. Our numerics break down at low temperatures due to our 
choice of gauge fixing, and we have stopped our numerics in a regime where the results 
are still reliable. Extrapolating further, it is possible that a singularity or instability arises 
at some finite temperature. However, what cannot happen is that we end up at a smooth 

3 Holding B and p fixed independently is not meaningful, as they are dimensionful parameters and are 
thus changed by a scale transformation. 
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Figure 3: Plot of the entropy versus temperature at fixed B 3 /p 2 = and B 3 /p 2 = .15 ±.002, 
for k = 2/\/3 (supersymmetric value). The numerical results show that a small B field 
causes a large drop in the entropy at low temperature. 



finite entropy extremal black hole, as we will show that no such solution exists at nonzero 
B (within our Ansatz, which assumes such properties as translation invariance). 

In the k < 1 and k = 1 cases the situation is dramatically different, as the entropy 
appears to go to a finite value (which depends on the value of B 3 /p 2 ), as shown for k = in 
Fig. 4. In the k — 1 case the thermodynamics is governed by spacelike warped AdS3 black 
hole solutions, which represent exact near-horizon geometries for our theory. At strictly zero 
temperature, for both k < 1 and k — 1, the full interpolating solutions acquire a relatively 
mild singularity at the horizon, unless B = in which case we recover the Reissner-Nordstrom 
solution. These singularities can be understood from a perturbative analysis. Thus, strictly 
speaking, in these cases there is no smooth extremal finite entropy geometry, just as there 
was not in the k > 1 case. However, it may be more physically relevant to focus on the 
behavior for small but nonzero temperature, in which case a residual entropy is evident. 

Besides our numerical results, which are valid for arbitrary B and p, we have carried out 
a perturbative analysis of the solutions valid to order B 2 . These solutions can be obtained 
analytically by methods analogous to those employed in the AdS/fluid dynamics literature, 
starting with [32]. Insofar as they overlap, the perturbative results corroborate our numerical 
findings. 

The organization of the remainder of this paper is as follows. In section 2, we briefly 
review the D = 5 Einstein/Maxwell theory with Chern-Simons term, including the definition 
of the boundary current and stress tensor for asymptotically AdSs solutions. In section 3 we 
present the Ansatz for uniform electric charge density and constant magnetic field, and derive 
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Figure 4: Plot of the entropy versus temperature at fixed B 3 /p 2 = and B 3 /p 2 = 10.5 ± .6, 
for k = 0. The entropy appears to go to a finite value (but see the text for comments on the 
strict zero temperature limit). 



the reduced field equations. In section 4, regularity and boundary conditions are discussed 
both at the horizon and in the asymptotic AdS 5 region. Standard AdS/CFT formalism is 
used to express physical quantities such as the entropy, chemical potential, currents and 
stress tensor in terms of these asymptotic data. In section 5, the near-horizon geometries 
of our solutions are constructed analytically, and we discuss the existence of solutions that 
interpolate between these and AdSs. In section 6, a perturbative expansion in powers of B 
is shown to be smooth, except for extremal solutions. In section 7, numerical results are 
presented, and a picture of the phase diagram is assembled. In section 8 we conclude with 
a discussion. A number of appendices include details of computations omitted in the main 
text, a general discussion on the existence of factorized solutions, and the construction of 
solutions with a near-horizon extremal BTZ factor. 



2 Action, Current, and Stress Tensor 

The action for 5-dimensional Einstein-Maxwell theory with negative cosmological constant, 
and Chern-Simons term, is given by 4 

Sem = J d 5 x^ (r + F MN F MN - ^) + S cs + S hndy (2.1) 



Conventions: R\„ K = d K T^ - d v T* K + 1^1*, - I^r*, and R„ v = R 



8 



where the Chern-Simons action is given by, 5 



s <*=i^gJ aafaf (2 ' 2 > 

For the value k = k s = 2/y / 3, the action coincides with the bosonic part of D = 5 minimal 
gauged supergravity. In this paper, however, k will often be kept general, thus allowing for 
values different from k s as well. Boundary terms in the action are required for the proper 
renormalization of various physical quantities [33, 34, 35]. In a coordinate system (r, x^) 
where g rfJl = asymptotically for fx — 0,1, 2, 3, the boundary action Sbndy is given by, 

*-* = sk L d ' x ^( K - 1 + i«w + 1 ( ln z) F ^) < 2 - 3 > 

Here, 7^ is the metric induced by qmn on the boundary, and K is the trace (with respect to 
7) of the extrinsic curvature of the boundary given by = {d r r y tlv )/i2 y /g^). Henceforth 
we set the AdS radius to unity: L — 1. The non-diffeomorphism invariant lnr term in 
the boundary action is needed to remove the divergence associated with the trace anomaly 

The Bianchi identity is dF = 0, while the field equations are given by, 

= d* F + kF A F 
Rmn = ^9mn + ~^F PQ FpqQmn — ^FmpFn P (2.4) 

2.1 Boundary Current and Stress Tensor 

For large r the boundary metric of asymptotically AdSs solutions will behave as 

7^ = r 2 ^J +■■ (2.5) 

where 7^ is the conformal boundary metric, given here by the flat Minkowski metric. Sim- 
ilarly, the components of the gauge field tangent to the boundary will go to a constant 
at large r, representing a constant magnetic field pointing in the x 3 direction. 

By considering on-shell variations of the boundary metric and gauge field, we can define 
a boundary current and stress tensor in the familiar fashion: 

5S = J d'x^W (j^A, + V^ 7 i° } ) (2.6) 



5 Note that our convention for k differs from that in [23]: fch crc = f^thcrc- The present convention has 
been chosen to simplify the Maxwell equation. By sign-reversal of A, we are free to choose k > 0. 
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Specializing to the case of a constant field strength on a flat boundary metric, we have 

-4irG 5 r = (r^F ri/ + ^e a ^A a F 01 ) (2.7) 

8nG 5 T^ = r^-K^ + K^ v - 3 7 ^ - 2 (F» a F u a - ^F a/3 F Q/37 ^ In r 

where the limit of large r is implied. For our solutions, the explicit lnr terms will cancel 
logarithmic terms in the metric functions, yielding a finite large r limit for T^ v . 



3 Ansatz and Reduced Field Equations 

The presence of uniform magnetic field and electric charge density in the boundary CFT 
may be achieved by an Ansatz for the bulk fields which is invariant under translations in x^, 
and space-rotations around the magnetic field, which we shall take to be pointing in the x 3 
direction. The Ansatz consistent with these symmetries is given by 6 

F = E(r)dr Adt + Bdx 1 A dx 2 + P(r)dx 3 A dr (3.1) 

for the Maxwell field strength, and by 

ds 2 = ^- U(r)dt 2 + e 2V » (dx\ + dxj) + e 2W ^ (dx 3 + C{r)dtf (3.2) 

for the metric. The magnetic field B is forced to be constant by the Bianchi identity, and 
the functions E, P, U, V, W, and C depend only on r. Reparametrization invariance in r has 
been used to select a coordinate r for which the same function U (r) appears in both g rr and 
the first factor of dt 2 . Rescaling x 1: x 2 can be compensated by a constant shift of V, while 
rescaling x 3 can be compensated by scaling C and shifting W. 

In this coordinate system the event horizon is located at U(r+) = 0, and the Hawking 
temperature is given by 

(3.3) 

An 

The above Ansatz is covariant under the following transformation, 

X3 x 3 — at 
C{r) -> C{r) + a 

E(r) -> E{r)-aP{r) (3.4) 



6 An additional term of the form N(r)dx 3 A dt must have constant N in view of the Bianchi identities, 
and N = in view of the field equations, and has thus been omitted. 
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with all other coordinates and fields left unchanged, and for any real parameter a. We shall 
refer to this transformation as asymmetry. It may be used, for example, to set C(oo) to 
zero. Note that the combinations E + CP and C are invariant under a-symmetry. 

The Ansatz (3.1), (3.2) is also covariant under boosts in the direction of the magnetic 
field. We can use these boosts to put the solution in the rest frame, defined by C{r + ) = 0. 
In the dual CFT this corresponds to setting to zero the chemical potential conjugate to 
momentum. See appendix A for the details. 

3.1 The reduced Maxwell equations 

The Bianchi identity is automatic on the Ansatz (3.1). The reduced Maxwell equations are, 

((E + CP)e 2V+w )' + 2kBP = 
(yPe 2V ~ w -C(E + CP)e 2V+w )' + 2kBE = (3.5) 

Both equations may be recast in terms of the a-symmetry invariant combinations C and 

S = E + PC (3.6) 

only, by eliminating the derivative of (E + CP)e 2V+w in the second equation using the first 
equation. One then obtains an alternative form of the reduced Maxwell equations, 

Ml (Se 2V+w )' ' + 2kBP = 

M2 (UPe 2V ~ w )' - C'£e 2V+w + 2kBS = (3.7) 

3.2 The reduced Einstein equations 

The reduced Einstein equations may be simplified to assume the following final form, 



El 



(C> e 2V+3 W y = A p Se 2V + W (38) 



E2 U(V" - W") + (U' + U(2V + W')) (V - W) 

= \ C 'fe 2W - 2B 2 e~ 4V + 2UP 2 e~ 2W 

E3 UV" + U'V + UV'(2V + W') = A- -S 2 - -B 2 e~ AV + -UP 2 e~ 2W 

3 3 3 

EA U" + U'(2V + W) - (C') 2 e 2W = 8+-£ 2 + -B 2 e^ v + -UP 2 e' 2W 

3 3 3 
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along with the constraint equation, 

CON U'(2V' + W) + 2U(V') 2 + AUV'W + ^(C') 2 e 2W 

= 12- 2S 2 - 2B 2 e~ 4V + 2UP 2 e - 2W (3.9) 

The r-derivative of the constraint vanishes by the use of the other six equations, and may 
be enforced as an initial condition, as usual. 

4 Asymptotics and initial data 

The solutions we consider are asymptotically AdS 5 . Thus, U(r), e 2V ^ and e 2W ^ behave as 
r 2 in the limit r — > oo. The precise overall normalization depends on the normalization of the 
space coordinates xi,x 2 (for V) and £3 (for W). There are two natural ways of normalizing 
this behavior, by parametrizing either the initial data at the horizon, or the asymptotic 
behavior as r — > 00. We shall consider both, and relate their behaviors. The data at the 
horizon is used for the numerical analysis, and for computing the entropy and temperature. 
The asymptotics are used to calculate the Maxwell current and the stress tensor. 



4.1 Parametrizing the initial data at the horizon 

In the numerical analysis it is important to choose coordinates to remove the gauge freedom. 
This can be done by demanding that the solution take a canonical form at the horizon. By 
rescaling of Xi,x 2 ,Xs, and combining an a-transformation and a boost in the ^-direction, 
the field strength Fh and the metric ds 2 H at the horizon may be arranged to take the form, 

Fh = qdr Adt + b dx\ A dx2 

ds 2 H = dx\ + dx\ + dx\ (4.1) 

where q and b are respectively the charge density and the magnetic field at the horizon (in the 
coordinates xi, x 2 , x 3 ). This corresponds to the following initial conditions at the horizon, 

£(r + ) = q U{r + ) = V(r + ) = W(r + ) = C{r + ) = P{r + ) = (4.2) 

We will refer to these coordinates as the horizon frame. It remains to specify V'(r + ), W'(r + ), 
and C'{r- + ). These quantities are generally not independent, but rather follow from the 
reduced equations M2, E2, and CON evaluated at the horizon, and we have, 

q (C"(r+) - 2kb) = 

U'(r + )V'(r + ) = A- 2 -q 2 - A -b 2 

U'(r + )W'(r + ) = ±-\{q 2 -b 2 )- l -C'{r + ) 2 (4.3) 
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The value of C'(r+) is specified to be C'(r + ) = 2kb for q ^ 0, but remains an independent 
free parameter for q = 0. The quantity U'(r + ) = 4ttT is not a genuine initial datum, since 
the equation CON for U is of first order. Therefore, genuinely distinct solutions are specified 
by only two parameters, for example T and q in units of magnetic field b. 

If the temperature T is nonzero, we can always rescale t to set U'(r + ) = 1, leaving the 
free parameters b and q. Furthermore, we can shift r to set the horizon at r + = 1. 



4.2 Extremal solutions require bq(k ± 1) = 

It is now easy to establish a non-existence result that will play an important role in what 
follows. We ask under what conditions can we have an extremal horizon, U (r+) = U'(r+) = 0. 

Assuming that all functions in our Ansatz are well behaved at r+, we can always work 
in the horizon frame specified in the previous subsection, in which case the conditions (4.3) 
apply. But then it is easy to see that the assumption of an extremal horizon with nonzero 
b and q is inconsistent with (4.3) unless k 2 = 1. To obtain an extremal horizon we must 
choose one (or more) of q — 0, b — 0, or k — ±1. As will be discussed in section 5, these 
three choices lead to near-horizon geometries of the form AdS2 x R 3 , AdS^ x R 2 , and warped 
AdS 3 x R 2 . But in the generic case in which none of these conditions is satisfied, finite 
temperature solutions cannot be smoothly brought to zero temperature, a feature that we 
will see explicitly from various points of view. 



4.3 Asymptotic behavior as r — > oo 

Starting from the initial data (4.2) and integrating our to large r, we will find asymptotically 
AdS 5 solutions with large r behavior 



U2 ,lnr _ c 4 

U=(r-r [ f + ^ +<— + ■■■ C = c + -^ 



e =v(r-r ) + - + v 2 — + ■ ■ ■ £ = ^ + 

where the dots stand for higher order terms in 1/r. Some of the parameters are related to 
one another by the field equations, 

2wv 2 , 2b 2 

W2 = u 2 = - — 

V 6V Z 

V2= 3v W2 = -^ (45) 

In these coordinates the conformal boundary metric is — dt 2 +v (dxl+dxl)+w(dx3+codt) 2 . As 
we show in appendix B, we can always perform a coordinate transformation to set v — w — 1 
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and Co = 0, which brings the conformal boundary metric to the standard Minkowski form. 
Further, this can be done while preserving the condition that we be in the rest frame, 
defined by C(r + ) = 0. The components of the current in this frame, which we refer to as 
the asymptotic frame since it is the relevant one for comparing with the boundary CFT, are 
computed to be 

2kb 

AnGsJ* = p = 7 c (e 3 - cop 3 ) - — (4.6) 

Sv 

J 1 ' 2 = 



47rG 5 J 3 = ( -^L - y/wc e 3 ) 

4 yy/w J 



Here B is the value of the magnetic field in the asymptotic frame, thus identified as the 
magnetic field in the CFT, and given by 

B = - (4.7) 

v 

7 C is a Lorentz boost factor (we use the notation 7 C to avoid confusion with the boundary 
metric 7), appearing when we transform to the rest frame, 

7c = (4.8) 

' 1 — WCn 



For black hole solutions in which g 33 remains finite at the horizon, A^l^ is arbitrary, 
and can be adjusted by a constant shift of A 3 throughout spacetime. Its appearance in 
the expression for J* is a consequence of the anomaly equation for the boundary current. 7 
A similar factor of A t \oo appears implicitly in J 3 , but since we have fixed A t = at the 
horizon, the asymptotic value of A t is determined without ambiguity. A nonzero value of 
A 3 |oo corresponds in the CFT to adding a chemical potential for J 3 ; it is simplest to set it 
to zero, and we do that henceforth unless stated otherwise. 

We similarly have expressions for the temperature, entropy density, and chemical poten- 
tial in the asymptotic frame: 

T = lcU ' {T+) (4.9) 

47T 



= — 



3lcV ( !- P3 " 



7 Given A 3 (t) and a constant magnetic field along x 3 we have d t J* <~ kE ■ B ~ kBd t A 3 . Integrating gives 
J* - kBA 3 , in accord with (2.7). 
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Comparing (4.6) and (4.9) we note the simple relation 

47rG 5 J 3 = -\kBn (4.10) 

This is the chiral magnetic current [13, 14, 15]: in the presence of a nonzero anomaly coeffi- 
cient k and a chemical potential, a current is induced parallel to the applied magnetic field. 
This result also follows from the anomaly equation, as can be seen by allowing for a slow 
variation in x^, differentiating both sides, and identifying V/i with an electric field. 

4.3.1 Physical quantities 

The global AdS$ solution is invariant under scale transformations, 

x'^ = x fM /£, r' = £r (4.11) 

Asymptotically AdS§ solutions inherit this as an asymptotic symmetry, reflecting the CFT 
nature of the holographic dual theory. Individual quantities, such as B and p, transform 
under these scalings, just as the coordinates x^ do, and so have no independent meaning. 
We should instead look at scale invariant quantities, which have physical meaning. This is 
the same as looking at dimensionless quantities from the boundary point of view. Let's write 
O ~ £ p to denote the transformation O' = £ p O. From the asymptotic invariance of the field 
strength F and the metric ds 2 , we find, 

s~£ 3 , T~£, B~£ 2 , P~£ 3 , p~£ (4.12) 

Any combination behaving as £° is a good physical quantity to compute. 

5 Near-horizon geometries 

For generic assignments of the physical parameters p and B, analytical solutions are not 
available in AdS$ (in contrast with the AdS^ case where an electric-magnetic duality rotation 
acting on the B = solution produces a simple dyonic solution). Even the special case of 
p = at zero temperature (and 5^0) does not lend itself to a full analytical solution [23]. 

Considerable qualitative and quantitative progress can be made, however, by solving for 
the near-horizon geometry of the solutions. This will be carried out in this section. Especially 
important will be the question as to whether, for given p, B, solutions with extremal near- 
horizon geometry exist, and whether they can support an electric field at the horizon. The 
existence of these extremal solutions is key to understanding the low temperature limit. One 
important result, which was already established in section 4.2 will be that for nonzero p and 
B, and for k ^ 1 there do not exist smooth, finite entropy, extremal solutions. 
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In the low temperature regime, the full solutions may then then be viewed as interpola- 
tions between asymptotic AdS§ and these near-horizon geometries. By numerical study, to 
be discussed in full in section 7, the regularity of this interpolation will be verified (except 
at strictly zero temperature, where singularities develop), and the physical properties of the 
solution, such as entropy, temperature, and mass will be evaluated. 

A general discussion on the existence of factorized solutions may be found in appendix 

C. 

5.1 General conditions for the existence of extremal solutions 

At an extremal horizon r + we have U(r + ) = U'(r + ) = 0. Extremal solutions provide a 
natural boundary of the parameter space of all solutions. 

To study their existence systematically, it will be convenient to adopt the horizon frame 
specified in section 4.1. We scale the coordinates Xj so that V(r + ) = W(r + ) = 0, and denote 
the magnetic field in these coordinates by b. Reduced Einstein/Maxwell equations M2, E2, 
E3, and CON, in which only U and U' enter on the left hand side, produce a set of non-trivial 
constraints, 

M2 q(C'(r + ) -2kb) = 

E2 C\r + f - Ab 2 = 

E3 6-q 2 -2b 2 = (5.1) 

where q is defined to be the electric field at the horizon q = £(r+). The constraint equation 
CON is a consequence of E2 and E3, and thus has been omitted from the above list. Elimi- 
nating C'(r + ) using the second equation reduces the system to qb(k±l) = and q 2 + 2b 2 = 6. 
The solutions are as follows, 

1. If k 7^ ±1, then we have qb = 0, so that either q = or b = 0, leading to the solutions, 

(a) The case 6 = and q = ±\/6, corresponds to the well-known extremal electrically 
charged black brane (without magnetic field, and arbitrary value of k). Its near- 
horizon geometry is AdS 2 x R 3 . 

(b) The case q = and b = ±\/3, corresponds to the extremal purely magnetic 
brane (without electric charge, and arbitrary value of k), obtained in [23]. Its 
near- horizon geometry is AdS^ x R 2 . 

2. If k — ±1, we shall show below that there is in fact a regular solution for every 
assignment satisfying q 2 + 2b 2 = 6, whose near-horizon geometry smoothly interpolates 
between AdS 2 x R 3 (at 6 = 0) and AdS 3 x R 2 (at q — 0). These solutions can be 
generalized to include finite temperature and momentum, and correspond precisely to 
one family of warped black holes considered in [30] . 
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5.2 Vanishing magnetic field: AdS2 x R 3 

We begin by briefly reviewing the well-known black brane solution in AdS 5 for B = 0, charge 
density p and mass M > (the actual charge and mass densities are proportional to these), 
given by the functions P = C = 0, and 

8 = ^ V = W = lnr U = r 2 + £ - ^ (5.2) 

In terms of the radii r_ < r + of the inner and outer horizons, we obtain a convenient 
parametrization of p, M and U, 

p 2 = 3r 2 r!(r 2 + r 2 ) 
M = rl + ri+r 2 + r 2 _ 

U = ^(r 2 - r 2 + )(r 2 - r 2 _)(r 2 + r\ + r 2 _) (5.3) 

As long as M 3 /p 4 > 3/4, the singularity at r = is protected by a horizon. The near-horizon 
metric reduces to 

dv 2 

ds 2 H = — — —— - U H {r-)dt 2 + r 2 + (dx\ + dx\ + dx\) (5.4) 
U H {r) 

with U H (r) ~ (r 2 — r 2 ) for the non-extremal case, and U H (r) ~ (r 2 — r^) 2 for the extremal 
case. The near-horizon geometry is factorized into the space-part which is flat, and an AdS2 
factor. The temperature T and entropy density s of the black brane are given by 

= (rj-r 2 _)(2rl + r 2 ) r^. 

2vrr3 4 1 ' J 

The black brane becomes extremal as r_ — > r + , so that the temperature goes to zero, but 
the charge density and entropy density remain finite, and related by s/p = l/(4\/6). 

5.3 Vanishing charge density: AdS% x R 2 

With vanishing charge density, the Maxwell field strength F reduces to the B-term only. 
The solutions in this case were obtained in [23]. The Maxwell-Einstein equations have an 
analytical solution, given by £ = P = C = 0, and 

f/ = 3(r 2 -r 2 ) e2V = ^ e 2W = r 2 (5.6) 

which represents the product of a (non-rotating) BTZ black hole with R 2 . It was confirmed 
numerically in [23] that there exists a family of regular solutions, parametrized by T / \J~B, 
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which interpolate between the BTZ black hole of (5.6) at the horizon, and AdS& at r = oo. 
The entropy of these solutions tends to zero as T — > 0, while the physical magnetic field B 
is kept fixed. 

More generally, we can have nonextremal rotating BTZ black holes, whose metric is given 
by e 2V = B/y/3, and 

U = 3 {r2 - r + ){ : 2 e™ = r 2 C = -V3 r -^r (5.7) 

with £ = P = 0. A useful alternative parametrization of the rotating BTZ solution is given 
by 

U = 12(r-r+)(r-r_) C = 2V3(r-r+) (5.8) 

and V = W = £ = P = 0. Note that both forms admit a smooth extremal limit. 



5.4 k = 1: warped AdS3 black holes 

As shown above, there is a special value of the Chern-Simons coupling, namely k = ±1 
(recall that this is less than the value required for supersymmetry), for which there exist 
extremal solutions for any q, b satisfying the relation, 

q 2 + 2b 2 = 6 (5.9) 

Furthermore, there is a simple nonextremal generalization. The solution is given by V = 
W = P = 0, and 

8 = q (5.10) 
U = 12(r — r + )(r — r_) 

C = 2b(r-r + ) (5.11) 

We have used a-symmetry to set C(r+) = 0. The metric and field strength are then given 
by 

dv^ / \ 2 

ds 2 = — — -U(r)dt 2 + (dx 3 + 2b(r - r+)dt) +dx\ + dx 2 2 (5.12) 
U[r) v ' 

F = qdr A dt + b dx 1 A dx 2 

The extremal limit is given by taking r_ = r + as usual. 

These solutions can be identified with the "self-dual" solutions described in section 6.1.1 
of [30], where we make the identification 

9 3b 2 , 
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Note that as b 2 ranges over its allowed values between and 3, v 2 ranges between and 
1. The equivalence can be seen most easily by comparing our metric with eqtn. 1.2 of [36]. 
Under the identifications 

12/ 2 / 12/2 

the metrics are seen to be proportional. Unlike in [36], we do not compactify $ ~ t. Our 
metric has no closed timelike curves or other pathologies. 

These solutions also arise in a context closely related to ours, namely 2 + 1 dimensional 
Einstein-Maxwell-Chern-Simons theory (see [37] for solutions which are the analytic contin- 
uation of these). This can be understood from the fact that if we reduce our theory down 
to three dimensions along x± i2 , then we recover the equations of 2 + 1 dimensional Einstein- 
Maxwell-Chern-Simons theory coupled to a massless scalar field. The condition for the scalar 
field to take a constant value, representing a constant value of V, is precisely the condition 
q 2 + 2b 2 = 6 that we found above. It is curious that these solutions exist only at the special 
value k = ±1. 

5.5 Existence of interpolating solutions 

As we have seen, for nonzero values of p and B, and k ^ ±1, there do not exist smooth 
zero temperature solutions under the assumptions of our Ansatz (it is possible that such 
solutions do exist if one, for example, relaxes the condition of translation invariance). So in 
these cases, if we start from an asymptotically AdS 5 solution at finite temperature, as we 
lower the temperature some of the functions in our solution will start to diverge; we will see 
this breakdown of our numerics. 

This leaves the question of what happens in the zero temperature limit in the cases for 
which there do exist candidate extremal horizons. In the case of B — and nonzero p 
(the value of k is immaterial in this case), the answer is that we end up at the usual AdS 5 
extremal Reissner-Nordstrom solution. This will turn out to be the only case in which we 
find a truly nonsingular extremal solution. 

For p = and nonzero B (k again drops out of the discussion if one notes that for q = 
the equation M2 appearing in (5.1) becomes trivial ) interpolating solutions were constructed 
numerically in [23]. At low temperatures the near-horizon geometry approaches AdS3 x R 2 , 
but as will be discussed momentarily a singularity develops in the full interpolating solution 
at strictly zero temperature. 

For nonzero p and B, but k = ±1, we have candidate near-horizon geometries corre- 
sponding to warped AdS3 black holes times R 2 , and these have a smooth extremal limit. 
At any nonzero temperature, our numerics will establish the existence of solutions smoothly 
interpolating between these near-horizon geometries and AdS 5 . As the temperature is taken 
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to zero the entropy remains finite, but nevertheless a singularity develops at the horizon, for 
reasons that can be seen as follows. 

To construct a candidate extremal interpolating solution, we can start with the exact near 
horizon extremal warped black hole geometry, and then introduce a perturbation that grows 
near the boundary, representing the change in asymptotic boundary conditions taking us 
towards AdSs. This perturbation is obtained by solving the equations obtained by linearizing 
around the near-horizon solution. 

For our near-horizon solution we have 

U = 12(r-l) 2 (5.15) 

C = 26(r-l) 

S = q 

V = W=P=0 



with q 2 + 2b 2 = 6. Note that we have used the freedom to rescale r to set the horizon at 
r + = 1. 

Now we linearize around this solution, denoting the perturbations by lower case letters (e 
denotes the perturbation of £). Plugging in we find the following equations to linear order 

Ml : 5e + q(2v' + w') + 2bp = (5. 16) 

M2: (l2(r- lfp) - qd - 2qb(2v + w) = 

El : c" + 2b(3w' + 2v') - Aqp = 

E2 : (l2(r - 1) V - w'j)' - 8b 2 v - 2bc - Ab 2 w = 

E3: (l2(r - 1)V)' - % 2 v + = 

EA : 3u" + 72(r - l)(2v' - w') + %b 2 {2v - 3w) - 126c - 16qe = 
CON : 24(r - l)(2v' + w') + Ab 2 (w - 2v) + 26c + Aqe = 

We impose the boundary conditions 

u(l) = u\l) = v(l) = w(l) = c(l) = c'(l) = p(l) = e(l) = (5.17) 

It is fairly straightforward to solve these equations iteratively. Of most relevance are the 
resulting expressions for v, e and p, which are 

v = 0l (r-l) a (5.18) 
e = 3qa 1 (r-l) a 

3 ?aai ^ q 

P = "^r (r - 1} "26 a2 
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where are integration constants, and 



1 V81 - 8b 2 

If q 2 7^ 6 then b 2 > 0, in which case a < 1, and then we see that p diverges at the horizon. 
Thus to have a smooth solution we are forced to set a\ = 0. But this means that v — e — 0, 
and from here it follows rapidly that the entire solution is just the original near-horizon 
geometry we started from. 

That the divergence in p represents a physical singularity can be seen by considering an 
infalling observer. One finds that such an observer sees a diverging physical field strength 
at the horizon. 

We therefore do not expect to find a smooth extremal interpolating solution when q ^ 0, 
even for k — ±1 when smooth near-horizon geometries do exist. In our numerics, as we 
lower the temperature we indeed find that P begins to diverge at the horizon in precisely 
the manner described above. Nevertheless, since the metric components have well defined 
limits (though not their derivatives) we find that the entropy appears to smoothly approach 
a finite value. 

This leaves the case of nonzero B but p = 0, which was studied in [23]. There solutions 
were found numerically that interpolated between near-horizon BTZ xR 2 and AdSs. At any 
finite temperature these interpolating solutions are smooth, but it follows from the above 
that a mild singularity develops at strictly zero temperature. For these solutions k does 
not appear in the field equations, and so we can freely set k — ±1, in which case we can 
compare with the linearized analysis just described by setting q = 0. For q = we have that 
p — 0, and so we avoid the divergence in that quantity. However, we have that v ~ (r — l) a , 
a = —1/2 + y/57/6. Since a < 1, the first derivative of V will diverge at the horizon, 
presumably indicating a singularity. 

Actually, one special case remains. When p = 0, rather than considering the zero temper- 
ature limit of BTZ corresponding to pure AdS3, we can look for solutions with near-horizon 
geometry given by a finite entropy extremal BTZ solution. These can be constructed by a 
fairly standard construction, as we describe in appendix D. The resulting interpolating solu- 
tions exhibit the same singularity at extremality as above; namely, derivatives of V diverge 
at the horizon. 



6 Perturbation theory in powers of B 

In this section, we shall construct solutions perturbatively in powers of B around the analyt- 
ically known solution for B = and arbitrary charge density p and mass M. We expect this 
expansion to be reliable for T 2 ^> B. For small T, however, the B = solution is close to 
extremal, and we know from the conditions of extremality of (5.1) that no extremal solutions 
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exist with B ^ and p ^ 0, unless we also have k = ±1. Thus, for k ^ ±1, we expect 
perturbation theory in B to break down near B ~ T 2 . For /c = ±1 the behavior is better, 
but we find that a mild singularity still results. The structure of the computation is very 
similar to the long wavelength fluid dynamics from gravity analysis in [32]. This is because 
a weak magnetic field corresponds to a slowly varying gauge field. 

The functions E(r), U(r), V(r), W(r) are even in B, while the functions P(r) and C(r) 
are odd in B. Here, we shall expand up to order B 2 included, so that, 

U = U + B 2 U 2 E = E + B 2 E 2 (6.1) 

V = V + B 2 V 2 C = Bd 

W = W + B 2 W 2 P = BP 1 

The zeroth order solution coincides with (5.2), and is given by P = C = 0, and 

E = ^ V = W = \nr U = r 2 + ^ - - (6.2) 

The horizons of Uq will be denoted by r±. The boundary conditions include Ci(r + ) = 
Ci(oo) = 0, together with the requirement that V 2 and W 2 fall off faster than r -2 as r — > oo. 
To separate out the spin zero (scalar) and spin two (tensor) perturbative corrections, it will 
be convenient to introduce the combinations, 

S 2 = 2V 2 + W 2 T 2 = V 2 - W 2 (6.3) 

Equations M2 and El are odd in B, and thus have contributions only to first order in B, 
while all other equations, Ml, E2, E3, E4, and CON are even in B and thus have only second 
order contributions. We begin by solving M2 and El first. 

6.1 Spin one sector 

To order B, equation M2 is given by, 

(6.4) 

It may be readily integrated, and the integration constant is fixed by the boundary conditions 
at the horizon, U (r + ) = Ci(r + ) = 0. As a result, Pi is given in terms of C\ by, 

Pi(r) = -jfr-, (Ci(r) + - 2 - 4) (6-5) 
rU (r) \ r 2 r\ J 
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Note that the function P\ is automatically smooth at the horizon. To obtain Ci, we substitute 
the solution Pi into equation El. Using the special form of the function Uq, this equation 
may be recast as follows, 

In this form, the equation may be solved by two successive integrations, producing two 
integration constants. These constants are fixed uniquely by the requirements that Ci(r+) = 
Ci(oo) = 0, and we find, 

Thus, the functions Pi and C\ are uniquely determined by the boundary conditions. 

6.2 Spin two sector 

Equation E2 to order B 2 gives, 

(r 3 U T>)' = l -r\C[f + 2rU Q (P 1 ) 2 - 2 - (6.8) 

The terms on the right hand side are known from the solution in the spin one sector. This 
equation may be solved by two successive integrations. The two resulting integration con- 
stants may be fixed by demanding smoothness of T' 2 at the horizon, and the vanishing of T 2 
at infinity, so that T 2 is uniquely given by, 

6.3 Spin zero sector 

The functions E%, Ui and ^2 all correspond to scalar perturbations. The linear combination 
3 x E3 — CON gives an equation for S2 in terms of C\ and T 2 , 

U Q S'i + 2 U S' 2 + U T% + 3 U T> + = Jr 2 (^) 2 - \ (6.10) 

r r Z r 

Eliminating T 2 using (6.9), the remaining equation can be recast in the form, (r 2 S' 2 y = —2P 2 , 
which may be solved by two successive integrations. We need S2 to fall off faster than 1/r 2 
at infinity to preserve the boundary. The solution is therefore uniquely fixed to be, 

S 2 (r) = 2 I" dr' f 1 - - ±) Pl (ry (6.11) 
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Next, equation Ml determines E 2 , 

(r 3 E 2 )' + (pS 2 )' + (r 3 P 1 C 1 )' + 2kP 1 = (6.12) 
The integration constant can be reabsorbed into p, and so we have 

E 2 (r) = -^S 2 - P 1 C 1 - ™ [dr> P 1 (r') (6.13) 

Equation E4 determines U 2 . It may be expressed as (r 3 U 2 )' = X, with X given by, 

(6.14) 



X(r) = -r 3 U^ 2 + 1 f(E 2 + C 1 P l )+r 5 (C' i y + 4rUoP? - 4 



3 + 3r 



The solution which goes to zero at infinity is 



r //3 i 
oo T Jr + 



dr'X(r') 



0-3 

2r2 



(6.15) 



where a 3 is an integration constant. Finally, the constraint equation may be checked to hold 
at r — > oo to leading order. This will guarantee that it is obeyed throughout. 



6.4 Asymptotic behavior of the perturbative solution 

The full perturbative solution is now fixed. The free parameters are: B, p, M, where we're 
not counting a 3 since it can be absorbed into M, nor r + since it is a function of M and p. 
As a result, the asymptotic behavior of these functions can now be computed, and we find, 
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6.5 Regularity of the perturbative expansion 



For the perturbative expansion around the non-extremal black brane with r_ < r+, the 
functions C\ and Pi fall off fast as r — > oo, and are smooth at the outer horizon r+, as well 
as at all other values of r > r+. As a result, the integrals giving S^, T 2 , -E2 and C/ 2 are rapidly 
convergent, and define regular functions throughout. 

The perturbative expansion around the extremal black brane with r_ = r+, however, is 
not, generally, well-behaved. For r_ = r+, the function C\(r) is given analytically by 



Ci(r) 



fc ^o(r) 

3 r 2 r 2 



2r 2 



r 2 + 2H 



1 , 

+ - In 
3 



r 2 + 2ri 



(6.17) 



The functions Ci and C( are smooth throughout. The function P\ is given analytically by 



Pi(r) 



kp 
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r 2 + 2r 2 



3r 2 + 2r 2 1 
+ 3(r 2 + 2r 2 ) + 3 



^2 ^,2 



r 2 + 2rl 



(6.18) 



and exhibits a logarithmic singularity as r — > r+. This singularity is integrable in the formulas 
giving the perturbation functions 5 , 2,T 2 , U2 and _E 2 . As a result, the functions 52, U2, E 2 are 
smooth, while the function T 2 has a logarithmic singularity as r — >■ r+, given by 



r 2 (r) 



67 + 



ln(r — r + ) 



(6.19) 



To summarize, in the extremal limit, the electric current density P, and the tensor pertur- 
bation of the metric T 2 both diverge logarithmically at the horizon, and the solution is not 
globally smooth. 



6.6 Perturbative calculation of entropy, temperature, and mass 

Having normalized the metric at r = 00 to be the conformally standard Minkowski metric, 
(B.l) with v — w — 1, the same metric at the horizon then reads, 

ds 2 H = r\ e 2V ^\di\ + dx\) + r 2 e 2W ^dx 2 3 (6.20) 

The perturbative corrections to the entropy density, temperature, and mass to order B 2 may 
then be deduced from the metric functions as follows. The entropy density s is given by 

s = s Q (l + B 2 S 2 (r + )) (6.21) 

Similarly, the temperature T and mass M are given by 

T = T + ^-B 2 U 2 (r + ) 

(12 lr°°/ 4 \ \ 

3-8 hr++ 2l*'( JfM -!sO) (6 22) 
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Here, sq, To, and Mq are respectively the entropy density, temperature and mass of the 
B — black brane. Note that it follows from the form of (6.11), that 5*2(r + ) < 0, so that 
the correction to the entropy density is always negative. 



6.6.1 Eliminating r + -dependence 

The position of the outer horizon, r + , has dimension and hence no direct physical meaning. 
It may be thought of as setting the overall scale, and plays the role of I above. Therefore, 
any physical quantity must be independent of r+. Thus, we shall introduce dimensionless 
coordinates and quantities, such as 

x = r/r + A = r_/r + (6.23) 

The r + -dependence may now be isolated in each one of the functions that enter perturbation 
theory. For C\ and Pi we define the dimensionless functions C\{x) and P\{x) by, 

d(r) = kp 2 C 1 {x)/r% Pi(r) = kpP 1 (x)/r 5 + (6.24) 

We also define the dimensionless functions er(A) and r(A) by, 

S 2 (r + ) = <j{\)/r\ U' 2 {r + ) = r(A)/rJ (6.25) 

and we use p as a physical quantity that sets the scale for r + , 

p 2 = r^(A) 6 z/(A) 6 = 3A 2 (1 + A 2 ) (6.26) 

so that r + = p 1 / 3 /z/(A). The expressions for entropy density and temperature normalized to 
the physical dimensionful quantity p are as follows, 

' / ;i-A 2 )(2 + A 2 ) + -^r(A)KA) 4 l (6.27) 



pV3 2ttz/(A) V 7V ' 2p 4 /3 

6.6.2 Calculating the dimensionless functions 

The dimensionless functions Ci(x) and Pi (a;) may be readily computed (analytically) from 
(6.7) and (6.5), and used to evaluate the dimensionless functions <r(A) and r(A), given by 

a(X) = -2k 2 v(\f J™ dx^-^P^x) 2 

r(A) = -^ + A; 2 A 2 (l + A 2 )f(A) (6.28) 
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where f (A) is given by 



f(A) = -^A 2 (l + A 2 )C^(1) 2 -9A 2 (1 + A 2 ) / dxx 3 C[(x) 2 -24 dxP^x) 

^i(z) + y cfeP^) 2 (8x - 4A 2 (1 + A 2 )^J (6.29) 

/oo „ /"OO / 1 \ ^ 

cfeP^x) 2 - 72A 2 (1 + A 2 ) J rfx M - -J P(:r) 2 

The functions er(A) and f(A) may be evaluated numerically, using the analytic expression for 
C\(x) and P\(x). The results are as follows. 

Numerical evaluation of the integral entering er(A) shows very little dependence on A 
throughout the interval A G [0, 1], and may be well approximated there by the average value 
of 0.015 (specifically, its value drops uniformly from 0.01505 at A = to 0.01490 at A = 1). 
As a result, we have the approximate formula, 

a{\) ~ -0.090 x A; 2 A 2 (1 + A 2 ) (6.30) 

for A throughout the interval [0, 1]. 

Numerical evaluation of the function f(A) produces a dependence given in figure 5 below. 
We record the end point values, 

f (0) = 8.605 r(0) = -1.333 

f(l) = 4.527 r(l) = -1.333 + 9.054 x k 2 (6.31) 

and for the supersymmetric value k 2 = 4/3, we have r(l) = 10.739. 

6.7 Physical interpretation of the perturbative corrections 

• If t(1) > 0, then T is non-zero and positive at the extremal value A = 1. We can extract 
an estimate for the minimum temperature under the assumption that B 3 /p 2 1, by 
simply estimating the temperature at A = 1, and we get, 

r{l)u{lf B 2 

Pmin = ~. (6.32) 

47T p 

which for the supersymmetric value k = 2/ \f?> gives approximately T min = 2.0932P 2 / p. 
Of course, at A = 1 our perturbative analysis breaks down, and so higher order terms 
could well invalidate this result. 



If r(l) < 0, then T = is attained for A = A c < 1, and the geometry must have a naked 
singularity whenever A c < A < 1, as it would correspond to negative temperature. 
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Figure 5: The function f (A). 

7 Numerical Analysis 
7.1 Setup 

We turn now to a discussion of our results obtained by numerical integration of the equations 
of motion. The first step is to specify our coordinate system. We impose the conditions 
corresponding to the horizon frame described in section (4.1), including choosing r + — 1 
and U'(l) = 1. This coordinate system will inevitably break down in the limit of vanishing 
temperature, since in that case we would have U'(r + ) = 0, and no rescaling of the time 
coordinate can bring us to our chosen gauge. We will see this breakdown occurring explicitly 
in the numerics. 

Solutions in this gauge are parameterized by the values of b and q, both of which we take 
to be non-negative without loss of generality. Choosing a value for the pair (b, q) fixes initial 
data at the horizon, and then we can integrate out to the asymptotically AdSs region at 
large r (we used Maple to do this). 

From the large r form of the obtained solution we can then compute the numerical 
coefficients (v, w, etc.) appearing in (4.4). We then convert these into physical quantities 
using the formulas given in section 4.3. The expressions we will be using in the following are 





b 


B = 






V 


T = 






An 




1 


s = 
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P = 7c(e3 - c p 3 ) 



with 7 C = 1/ — wcq. 

It is most illuminating to provide plots of entropy density versus temperature with the 
magnetic field and charge density held fixed. However, it only makes sense to keep fixed 
the dimensionless ratio B 3 /p 2 . Similarly, it is only meaningful to compute dimensionless 
versions of the entropy and temperature, and for these we choose 

(7.2) 



( p 2 + £3)1/2 ' (p2 + £3)1/6 

An instructive special case is the Reissner-Nordstrom solution with B = reviewed in 
section 5.2. The solution is originally given in terms of r + and p. Transforming this solution 
into the gauge used here, we find q = p/r+, along with 

8 -=± jL = L(^m (73) 

Note that the extremal limit in this parametrization is q = VE with s/p = l/(4\/6) « .102. 
A consistency check on our numerics is that we recover the curve described by (7.3) along 
with the extremal endpoint. 

The next step is to determine the region of the (b, q) parameter space that gives rise 
to smooth solutions. The boundary of this region depends on the value of k. Numerical 
integration shows that as we move out radially from the origin in the (b, q) plane we eventually 
find that some of the parameters v, w, etc. start to diverge or go to zero as we approach a 
curve in the (b, q) plane, depicted by the red lines in Fig. 2. The analytic form of this curve 
is only known at k — 1, where it is given by q 2 + 2b 2 = 6. Points on the k — 1 critical curve 
correspond to extremal warped AdS3 black holes, as was described in section (5.4). For the 
other values of k that we consider, the curve acquires some bulges, but continues to look 
roughly like that for k — 1. In the generic case, we determined the critical curve numerically 
by evolving outward until 7 C exceeds some specified value, which we take to be roughly 12. 
This requires fine tuning (b, q) to the critical curve to roughly four decimal places. The k = 
case is special since here 7 C = 1 exactly for all solutions; here we locate the critical curve by 
looking for a divergence in \ jv. In all cases, the dimensionless temperature tends to zero as 
we approach the critical curve. 

Once the critical curve has been identified, we can set up a grid in the (6, q) plane, and 
scan over the gridpoints. We took about 12, 000 roughly evenly spaced gridpoints. Once 
data at the gridpoints has been obtained we can search for curves along which B 3 /p 2 is 
approximately constant. These curves were illustrated schematically by the flow diagrams 
in Fig. 2. Finally, we plot the entropy versus temperature for points along such a curve. 
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7.2 Results 



We will discuss three cases: k — 0, k — 1, and k = 2/>/3, the latter being the supersymmetric 
value. We expect that these are representative of the general cases k < 0, k — 1 and fc > 1. 

Curves of approximately fixed B 3 /p 2 are shown in Fig. 6, which may be compared with 
the schematic version in Fig. 2. In the k = and k = 2/y/3 cases, the curves appear to 



k = 



k = 1 



10.5 ±.6 



w .1 ± .005 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.S 



0.2 0/4 0.6 0.8 



1.2 1.4 



A; = 



\ « .15 ± .002 
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Figure 6: Points in the (b, q) plane (b on horizontal axis) obtained by scanning over a grid 
and keeping points with B 3 / p 2 fixed within some interval. These numerical plots are to be 
compared with the schematic flows illustrated in Fig. 2 

be heading towards 6 = and q = respectively. As they do so, they begin to approach 
very near to the critical curve discussed above. We cannot follow them all the way there, as 
our limited precision prevents us from collecting data points arbitrarily close to the critical 
curves. It is conceivable that the true curves instead terminate at some location on the 
critical curves, as apparently occurs in the k — 1 case. 

Given the points along a fixed B 3 / p 2 curve, we can construct a plot of entropy versus 
temperature. Such plots for the k = 2/y/3 and k — cases are displayed in Figs. 3 and 4. 
In Fig. 7 we show the corresponding plot for k — 1. In all of these plots we compare the 
finite B results against those for B = 0. The B = curves represent the Reissner-Nordstrom 
black brane solution, and reproduce numerically the curve described in (7.3). 

In the k — 0, 1 cases, the entropy appears to go to a finite value at extremality, even 
though a singularity seems to be developing at the horizon in this limit, as was discussed 
in section 5.5. By adjusting B 3 /p 2 , we can tune this limiting entropy to any desired value 
between and that of the Reissner-Nordstrom solution at B — 0. In the plots we have chosen 
the value of B 3 / p 2 such that the extremal entropy is roughly half of its maximal value. We 
see that this requires a much larger magnetic field in the k = compared to k — 1. 

At k — 1 the entropy is controlled by the near-horizon warped AdS3 black hole solution. 

The situation for k > 1 appears to be very different, at least for the values that we have 
studied. The main effect, which becomes more pronounced at larger k, is that the entropy 
decreases substantially at low temperatures, and appears to be headed towards zero, until 
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S/{p 2 + B 3 )^ 2 




Figure 7: Entropy versus temperature for k — 1 and B 3 /p 2 .1 ± .005 



our numerics break down. For the supersymmetric case k = 2/y / 3 shown in Fig. 3 the effect 
is relatively modest due to the fact that 2/y/3 ~ 1.15 is not too much larger than 1. But 
even in this case it is evident that at low temperature a value of B 3 /p 2 much less than 1 
causes a decrease in the entropy by factor much larger than 1, and the effect grows as the 
temperature is decreased. 

The main property driving this behavior is the apparent location of the endpoint of the 
flows at fixed B 3 /p 2 . We can locate the endpoint numerically by setting q = and increasing 
b until a singularity (or numerical breakdown) occurs. Calling this endpoint value b c , we find 
b c = y/3 = 1.732 for k < 1, but b c < V% for k > 1. For k = 2/^3 we find b c « 1.568. 
On the other hand, with q = we showed that only for b = \/3 does a smooth, extremal, 
near-horizon geometry exist. Since there is no candidate smooth endpoint, it seems likely 
from this perspective that the flows terminate in a singularity at the point (q — 0, b — b c ). 
The metric functions we compute certainly behave very badly as this point is approached, 
but it is difficult to completely disentangle physical divergences from a breakdown of our 
numerics. 

The breakdown of our numerics at very low temperature is a reflection of our gauge fixing 
choice U'(r + ) = 1. Because of this, we are not able to unambiguously determine whether the 
k > 1 curves end at zero temperature and entropy, or terminate before then in a singularity. 
But we certainly do not see any evidence of a finite entropy endpoint at zero temperature, 
a conclusion which is bolstered both by our perturbative analysis and by the non-existence 
of any candidate near extremal geometries to match onto. The most likely scenario seems 
to be that finite entropy at extremality requires the magnetic field to be fine tuned to zero. 
One clear goal for the future is to improve the numerical treatment to allow a more refined 
study of the very low temperature regime. 
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8 Discussion 



In this work we have constructed asymptotically AdS 5 black brane solutions carrying nonzero 
charge density and magnetic field. They were found analytically in a perturbative expansion 
for small B, and numerically for general values. The most interesting results centered around 
the low temperature regime, where we found a sensitive dependence on both the magnetic 
field and the value of the Chern-Simons coupling k. For k > 1, including the supersymmetric 
value of k — 2/y / 3, we found that a small magnetic field causes a rapid decrease of the entropy 
at low temperatures. We proposed an AdS/CFT version of Nernst's theorem and the third 
law of thermodynamics consistent with the observed behavior. More general tests were left 
to the future [29]. 

k = 1 emerged as a special value, for here, and only here, there exist smooth finite en- 
tropy extremal near-horizon geometries carrying charge and magnetic field. We identified 
these, and their nonextremal generalizations, with one class of warped AdS 3 x R 2 black hole 
solutions studied in [30] (without the R 2 factor) in the context of topologically massive grav- 
ity. We found that these could be connected to asymptotically AdSs spacetimes, although a 
singularity in the interpolating solution develops at the horizon in the strict extremal limit. 

It is an intriguing question as to whether the value k = 1, which appears special from 
the point of view of supergravity solutions, has a special significance also on the CFT side, 
perhaps because it corresponds to a special embedding of the gauge field U(l) C SU(A) R . 

There are some other contexts in which sufficiently large values of the Chern-Simons 
coupling k causes novel effects. In the recent work [38] it was found that sufficiently large k 
can cause an instability in Reissner-Nordstrom black brane solutions. In [39] it was found 
that large k likely causes an instability of spinning black hole solutions. In both of these 
cases k must be larger than the supersymmetric value for an instability to occur, while we 
have seen here that the extremal Reissner-Nordstrom entropy is apparently destabilized even 
at the supersymmetric value. 

There should be many applications of these solutions to the study of condensed mat- 
ter and finite density QCD. Transport properties can be computed in these backgrounds; 
although we lack analytical solutions, a numerical treatment should be tractable. 

We conclude with a curious observation concerning the frequent appearance of 3/4 in this 
subject. To wit: the ratio of the high temperature entropy in gravity to that in the gauge 
theory is 3/4; the corresponding ratio of the low temperature entropy in the presence of a 
magnetic field is v/4/3; and the ratio of the "special" value of k (k — 1 in our conventions) 



to the supersymmetric value is y3/4. These three factors are not logically related in any 
obvious way. Perhaps this is just coincidence. 




32 



Note Added: 

Relaxing the gauge condition P(r + ) = by treating C'(r + ) as a tunable free parameter, 
and using numerical analysis with much higher precision than was used in the present paper, 
it has become possible to explore lower magnetic fields and lower temperatures in a reliable 
manner [42]. The results confirm that, for the supersymmetric value k = 2/y/3, the entropy 
density vanishes linearly with T when B 3 /p 2 > 0.124569, but reveals that for magnetic fields 
smaller than this value, the entropy density does not vanish at T = 0. It is established in 
[42] that the point T = 0, B 3 /p 2 ~ 0.124569 corresponds to a quantum critical point with 
dynamical scaling exponent z — 3. For any fixed B 3 /p 2 above the critical value 0.12459, 
and possibly for all values, the improved numerics also indicate that the flows towards zero 
temperature end at the AdS 3 x R 2 fixed point at b = y/3, q — 0, corresponding to the dot 
in Figure 2. It is expected that the low magnetic field phase should ultimately be unstable 
against turning on further couplings A«. 
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A Asymptotic boost symmetry 

The Ansatz (3.1), (3.2) is covariant under boosts in the direction of the magnetic field. 
Assuming the metric ds 2 to be asymptotically AdS 5 , the boundary space-time coordinates 
x M may be rescaled so that its asymptotic behavior, as r — > oo, is given by, 



produces a field strength and a metric of the same form as (3.1) and (3.2), but with the coor- 
dinates r, t, X3 replaced by f, t, x 3 and the functions E, P, U, V, W, C of r replaced by the func- 
tions E, P, U, V, W, C of f respectively, leaving x±, X2, and B unchanged. (A reparametriza- 
tion r — > f is generally needed to put the [/-function back into the gauge of the Ansatz.) 




(A.4) 



Performing a boost in the a^-direction 



%3 = lc(x 3 + /3t) 



7 C 2 (1-/? 2 ) = 1 
\0\ <1 



(A.5) 
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The relation between the transformed and original Maxwell fields are as follows, 

Edf = 7 C - /3P)dr 

Pdf = -f c (P-f3E)dr (A.6) 

while for the metric fields, we have \J~ l df = C/ _1 dr, and V = V, as well as, 

U-e 2W C 2 = - ^{C + pf e 2W 

e™ = ^(l + pcfe^-^U (A.7) 

and the following transformation law between C and C, 

.J C + 0)(l + 0C)^- fi U 
(l + f3C) e 2W - f3 2 U 

If C(oo) = 0, then it follows that C(oo) = for all (3. At an event horizon, r = r+, where 
U (r+) = 0, we have a simplified formula, 

As a result, in these coordinates, the value of C at the horizon characterizes the velocity of 
the configuration. Performing a boost j3 = — C(r+) brings the solution to its rest frame. 



B Relation between horizon and asymptotic frames 

Our solutions being asymptotically AdS$ implies that by rescaling t, x±, X2, £3, and perform- 
ing an a-transformation and a boost, the asymptotic behavior of the field strength and of 
the metric as r — > 00 may be put in standard form, 

F = Edr A dt + Bd5ti A dx 2 + Pdx 3 A dr 
ds 2 ~ ^ +r 2 (-dP + dxj + dxj + dxf) (B.l) 

In general, the coordinates t, X{ will be different from the coordinates t, Xi used in specifying 
initial conditions at the horizon. In the coordinates t,Xi, the asymptotics of the functions 
take the form displayed in (4.4) and (4.5). The rescaling and a-transformation relating the 
coordinates t, Xi and t, X{ is given by, 

t = t x lt2 = Vvx lj2 5t 3 = t/w (x 3 + cot) (B.2) 
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This is combined with a constant shift in the r coordinate by an amount tq. Since the 
magnetic field term in F is r-independent, we obtain a simple relation between the magnetic 
field b at the horizon, and the magnetic field B at infinity, 

b = vB (B.3) 

The functions U, V, W, C, £, P of the full metric in the coordinates t, Xi are then given by, 

2 u 2 2B 2 \nr ~ Jw c 4 

e 2 W+^ + ^ + ... £=5 + ... 

e «*=r>-*!|-?f^ + ... p = ^ + ... (B.4) 

While the metric at the boundary of ArfSs now takes the standard form of (B.l), the metric 
ds 2 H at the horizon has been rescaled and a-transformed to become, 

11 

ds 2 H = -{dx\ + dx\) H (dx 3 — ca\fwdt) 2 (B.5) 

v w 



corresponding to a solution moving with velocity Cq^/w. 

We can move to the rest frame of the solution by performing a boost in the ^-direction, 
thereby removing the cross term in ds 2 H while preserving the conformal boundary metric, 

t = 7 c (i + y/wcox?) , x 3 = 7 C (£ 3 + y/wcoi) , 7 C = (B.6) 

v 1 - wc l 

In these coordinates the field strength is 

- b ( P \ 

F = ^cEdf Adt H — dx.\ A c?a; 2 + —= \fwc^ c E dx 3 A df (B.7) 

w \ Vw7 c J 



and the metric on the horizon is 



ds 2 H = 1 (dx 1 2 + dx 2 ) + -^c£r 3 2 (B.8) 
v wy£ 



which corresponds to a solution at rest in the hatted frame. This hatted (asymptotic) frame 
is the one that we use to express physical quantities in the CFT. 
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B.l Evaluating entropy, temperature and chemical potential 



Now we can give expressions for the various physical quantities in these coordinates. The 
entropy density s is read off from metric on the horizon, and is given by 

s = \~^= 2 ( B - 9 ) 

The factors in the denominator ensure that s is the physical entropy density per unit volume, 
as measured in the CFT. 

To evaluate the temperature, we first recall its behavior under a boost. Let T be the 
temperature of some system in the rest frame. Then in a boosted frame the temperature is 
given by T = T /7 C , as can be seen by writing the Boltzmann factor e~ E ^ T ° in terms of the 
boosted energy and momentum. Now in the original (unhatted) frame the temperature is 
T = C/'(r+)/(47r). This represents the moving frame, while the rest frame is one in which 
C(r + ) = 0, which is the hatted frame. Therefore, the temperature in the rest frame, which 
is what we're interested in, is 

f = T = lc f= lcU '^ +) (B.10) 

From now on, when we write T we have in mind T. 

We also want the chemical potential, given by the asymptotic value of A t . Since we set 
A t = at the horizon we have A t \ 00 = J^df F rt . This gives 



fi = A = 7 C / 

oo J r 



dr E (B.ll) 

r+ 



If kB 7^ 0, we may use the second Maxwell equation, in the form given in (3.5) to recast E as 
a total r-derivative. The chemical potential may now be evaluated in the original unprimed 
coordinates, where the magnetic field is b, and we find, 



B.2 Evaluating the Maxwell current 

We identify the Maxwell current in terms of the on-shell variation of the Maxwell-Einstein 
action of (2.1), with respect to the Maxwell field at the boundary of AdS 5 , 

5S = Jd! L x\p)<®Ji i 8A li (B.13) 
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Here, 7^ is the asymptotic conformal metric defined by 7 Mi //T 2 — > 7^ in the limit r — > 00, 
where 7^ is the metric induced from g^ u on a surface of constant r. The on-shell variation 
of the action with respect to the gauge field may be calculated from (2.1) and we find the 
following expression for the current, 



4nG 5 J" = -(r 3 F r " + r 4 (\nr)d v F^ + ^e a ^A a F M ) (B.14) 

Here, e al3lfM denotes the volume form for 7^. Using the leading forms g rr = r 2 and g^ u = 
7(o) l r 2 we can wr ite this as 

4nG 5 J» = (r 3 ^F rv + (lnr)7^7(S^F a/9 + ^e a ^A a F M ) (B.15) 

Since we only consider solutions with constant field strength on the boundary the middle 
term will not contribute. 

B.3 Evaluating the current in the rest frame 

The current can be computed directly in the double primed coordinates system as 

4nG 5 J l = (e 3 -^=A 3 ) (B.16) 



3/ 



4nG 5 J ' = 

V w 3Vv 2 w 

The double primed coordinates have been defined so that v — w — 1. From the formulas 
given above we also have, 

B= b , e 3 = 7 c (e 3 - c p 3 ) , p 3 = 7 C ( - Vwc e 3 J (B.17) 

and 

We therefore finally have 



3 7^ ( r- \ 3u . , . 

Vwcoe 3 = -^ttP3 (B.18) 



2A;6 - 

47rG 5 J* = 7c(e 3 - c p 3 ) - ^3 (B.19) 

3t> 00 

47rG 5 J 3 = ^ c \^= - y/wc Q e 3 \ 
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to 



We note that A3 can be chosen arbitrarily as long as #33 is finite at the horizon (though we 
should set it to zero unless we want to add a source for J 3 in the CFT partition function), 

while A t is given in (B.ll). 

00 

C Factorized Solutions 

In the general 5-dimensional solutions we have obtained here, the 2-dimensional Xi,x 2 - 
plane perpendicular to the magnetic field is warped over the remaining 3-dimensional space. 
The near-horizon geometry of these solutions, however, invariably reduces to a space-time 
in which the xi, ^2-plane factorizes from the solutions because the field V becomes r- 
independent. This raises the question as to the structure of the most general factorized 
solution. In this appendix we shall show that, modulo certain regularity conditions, the only 
factorized solutions are of the type given in section 4. 

We begin by proving the following auxiliary result: if the function P vanishes, then we 
only have the factorized k 2 = 1 solution of section 4.4. Next, assuming now that P^O, but 
V = 0, we show that again only the factorized k 2 = 1 solution of section 4.4 exists. 



C.l P = leads to the near- horizon geometry 

By shifting V, W by constants, we may always assume that V(r + ) = W(r + ) = 0. Recall 
that the magnetic field in these coordinates is denoted by b. We assume kb 7^ and q 7^ 0. 
We begin by showing that the condition P = leads to U", V, W, £ constant. When P — 0, 
equations Ml, M2, El may be integrated, and give, 

£ = q e- 2V ~ w 
C = 2kbe- 2V ~ w 

C = 2kbe- 2V ~ 3W (C.l) 

Comparing the two expressions above for C, we find that W = identically. Taking the 
difference between E2 and E3 (for W = 0), we find, 

= 24-4£ 2 + Ab 2 e~ 4V -3(C) 2 (C.2) 

Substituting the above expression for C and £, we find, 

= 24 + e~ AV (-4g 2 + 46 2 - I2k 2 b 2 ) (C.3) 

Since the terms in the parenthesis are constants, V must be constant, and hence V = in 
view of the boundary condition. Thus, £ and C are constant. Equation E2 now requires 
= 2k 2 b 2 — 2b 2 which requires k 2 = 1. Equation E3 requires 

= 6 - q 2 - 2b 2 (C.4) 
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which is precisely the boundary curve equation for k 2 = 1. CON is automatic, and E4 gives 
U" = 24. But this gives precisely the factorized solution of section 4.4. 



C.2 V = leads to the near-horizon geometry 

Next, we assume P^O and V = 0, so that the geometry is factorized. Under this assump- 
tion, equation M2 may be traded for the constraint equation, since upon differentiation, CON 
will require the use of M2 when P ^ 0. Clearly, we must now retain E4 as an independent 
equation. Equations E2, E3, and CON are equivalent to the following equations, 



g 2 = 6 _ 2b 2 + UF 2 e -2W 

U'W = - l -{C'fe 2W + 2b 2 
-W"-(W') 2 = 2P 2 e- 2W (C.5) 

By eliminating S 2 , (C) 2 , and P 2 from E4, using the above equations, we obtain a relatively 
simple equation relating only U and W, 

U" + WW + 2UW" + 2U{W') 2 = 24 (C.6) 

Finally, equation Ml may be viewed as giving P in terms of £ and W, 

(£e w ) (C.7) 



P = 



2kb 

and El may be integrated upon using Ml, to give, 

C' = -±-£ 2 e- w + ae- 3W (C.8) 
kb 

where a is an integration constant. As promised, we omit equation M2. 

Equation (C.6) may be readily integrated in terms of the function r, defined by, 

~~dr'e w ^ (C.9) 



t r 



and we have 



U=^-rl) (CIO) 



where r and r are integration constants. Using the relations r' = e w and, 



t" t'" 



W' = — w" + (W') 2 = — (C.ll) 

T T 

39 



and eliminating £, U, P, W from (C.5) in favor of r and E = t'S, we find the following three 
equivalent equations, 



E 2 = (6-26 2 )(r') 2 -— (r 2 -r 2 ) (C.12) 

T 

(E'f = -2k 2 b 2 r'r'" 

= 48(r'0 2 (r 2 -r 2 )-48r(rOV + 46 2 (r0 4 -^(^ 2 -aA;6) 2 

Finally, E and E' may be eliminated from the above equations to yield two equations for 
t. Eliminating E 2 between the first and third equations yields an equation of third order 
in t. Eliminating E' requires some extra care, as we do not wish to unduly increase the 
order of the resulting equation, or introduce square roots. To obtain a second equation of 
third order, we take the product of the first two equations giving [(.E 2 )'] 2 /4, and use the last 
equation to eliminate E 2 from this. The results are as follows, 

= 6(6-2b 2 )r 2 (r"') 2 + 12r 2 b 2 (r"') 2 + (6-2b 2 ) 2 (r') 2 (r") 2 
-12(6 - 2b 2 )rr'r"r"' + 2b 2 (6 - 26 2 )(r')V" 



= A8k 2 b 2 (r") 2 (r 2 -r 2 )-A8k 2 b 2 T(r') 2 T" + Ak 2 b A (T') A 

- ^(6 - 26 2 )(r') 2 - 6^(r 2 - r 2 ) - q 2 - 2^ (C.13) 

We shall now examine the existence of joint solutions to both equations, which solve the 
boundary conditions, U(r+) = W(r + ) = P(r+) = 0, and C'(r + ) = 2kb. These conditions 
translate as follows 

r(r+) = r kba = q 2 + 2k 2 b 2 

r'(r + ) = 1 

r"'(r + ) = E'(r + )=0 (C.14) 

Evaluating the first equation of (C.12) at the horizon, using the above boundary conditions, 
gives the boundary curve relation, 

q 2 + 2B 2 = Q (C.15) 

Evaluating both equations of (C.13) at the horizon, and using the boundary conditions, 
r '(r + ) = 1 and r'"(r+) = 0, we find, 

= (6 - 26 2 )r // (r + ) 2 

= -A8k 2 b 2 r r"(r + )+Ak 2 b A -4:k A b 4 (C.16) 
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where in the last equation we have used the relation (C.15). Since we assume that q 7^ 0, 
relation (C.15) gives 6 — 2ft 2 7^ 0, and thus we must have r"(r + ) = from the first equation 
in (C.16). Using this result in the second equation gives altogether 

r"(r + ) = k 2 = l (C.17) 

Note that we have now, in principle, overdetermined the boundary conditions even just on 
a single equation, as we have r(r + ) = r ,r / (r + ) = l,r"(r + ) = r //; (r + ) = for a differential 
equation which is of third order. To proceed further, it appears necessary to make some 
assumption on the regularity properties of r near the horizon. A general Ansatz consistent 
with the above boundary conditions is as follows, 

r(r) = To + (r — r+) + r p (r — r + ) p + higher orders (C.18) 

for any real number p > 3. Substituting this Ansatz into the first equation of (C.13) shows 
that, as r — > r+, the last term is of order (r — r + )^ p ~ 3 ^ and dominates the other 4 terms, 
which all vanish faster as r — > r+. Thus, we must have r p = 0, for any p > 3. Thus, within 
the class of asymptotic behaviors given by (C.18), the expression 

r(r) = r + (r - r+) (C.19) 

is the only solution satisfying the boundary conditions. Clearly, for this solution, we have 
U" = 24, and W constant, which is the solution of section 4.4. 

D Interpolating between extremal BTZ xR 2 and AdSs 

In section 5.5 we alluded to the existence of solutions that interpolate between a near- 
horizon extremal BTZ xR 2 geometry and AdSs. In this appendix we give the details of 
their construction. They can be thought of as infinitely boosted versions of the solutions 
studied in [23]. 

D.l Extremal solutions with momentum 

In [23] we found zero temperature solutions interpolating between AdS3 x R 2 (with magnetic 
flux on the R 2 ) and AdSs. The metric and field strength are 

ds 2 = e-^d^ + e^i-d^ + dxD + e^idxl + dxj) (D.l) 
F = Bdxi A dx 2 

The Einstein-Maxwell equations reduce to 

2V" + W" + 2(V') 2 + {W') 2 = (D.2) 
(y>y + (w'f + 4 V 'W = 6e -w _ e -4V~2w B 2 
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One can (numerically) find solutions with small r behavior e 2V = _B/\/3, e 2W = r 2 , and large 
r behavior e 2V = vr 2 , e 2W = r 2 , with v ~ 1.87. Since B can be set to unity be a coordinate 
rescaling, there is actually a unique such solution. The solutions are not entirely smooth, as 
V develops a sub leading small r dependence r a with < a < 1, as was discussed in section 
5.5. 

Physically, these solutions represent the RG flow of M = 4 SYM theory in the presence 
of a magnetic field. At low energies the theory is governed by fermions in the lowest Landau 
level. These fermions are free to move parallel to the magnetic field lines, and give rise to a 
1 + 1 dimensional CFT at low energies, hence the appearance of AdS 3 . In [23] the central 
charges in gravity and M = 4 SYM at vanishing coupling were compared, and found to differ 
by a factor of y^3/4. Note that the theory is nonsupersymmetric even at zero temperature 
due to the presence of the magnetic field. 

On the CFT side, one should be able to excite one chirality of fermions to arrive at zero 
temperature configurations carrying momentum, and we expect corresponding solutions on 
the gravity side as well. The structure of such solutions follows, as in [40, 41], from the 
existence of null translational isometries. 

We take as our Ansatz 

ds 2 = e- 2W ^Ur 2 + e 2W ^\-dt 2 + dx 2 3 ) + e 2V ^(dx 2 1 + dx 2 2 )-u(r)(dx 3 -dt) 2 (D.3) 
F = Bdxi A dx2 

with V and W obeying (D.2). Plugging in, we find that the Einstein-Maxwell equations 
reduce to the following linear equation 

u" + (2V - W')u - 2{W" + {W'f + 2V'W')u = (D.4) 

This equation can be solved subject to the boundary condition that u(r) should fall off 
as 1/r 2 as r — > oo. The solution is given by 

u (r) = p e 2W ^ [ r di e - 2V ^- 3W ® (D.5) 

J oo 

where p is an integration constant. Given the asymptotics of V and W, the precise fall-off is 

«(r) = -£ + ••• (D.6) 
Evaluating the boundary stress tensor on these solutions we find 

8nG 5 T u = -\u 2 + l -p (D.7) 
87rG 5 T t3 = \p 
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where u<i determines the large r falloff of U — e 2W according to U = r 2 + + • • •. The form 
of the stress tensor is consistent with a momentum p being carried by chiral excitations. 
Another manifestation of this is the formula for the entropy density, which is 



1 „ fc 



Gs s = 2^-L (D.8) 
with 

The central charge is the same as was computed in [23] . 

The solutions just described carry momentum but no electric charge or current. The 
latter can be included in a fairly trivial way by turning on constant gauge potentials subject 
to regularity conditions at the horizon. According to (2.7), the combination of the Chern- 
Simons coupling k and the magnetic field B converts such constant potentials into charges 
and currents. The regularity condition at the horizon is A t + A 3 , and so we can turn on 
A t = —A 3 = A for any constant A. This induces 

47rG 5 J* = 4nG 5 J 3 = \kBA (D.10) 

Although we have thus constructed finite entropy, extremal, solutions carrying both 
charge and magnetic field, we see that this requires that we have a nonzero value of A 3 
at infinity, which corresponds to a nonzero chemical potential for current in the CFT. If 
we demand that this chemical potential is zero, corresponding to the boundary condition 
^3 |oo = 0, then these solutions do not appear. 
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